Spin dynamics of a one-dimensional spin-| fully anisotropic Ising-like 
antiferromagnet in a transverse magnetic field 



Asimkumar Ghosh 
Department of Physics, Scottish Church College, 1 & 3 Urquhart Square, 

Kolkata 700 006, India. 

Abstract 

We consider the one-dimensional Ising-like fully anisotropic S=^ Heisenberg 
antiferromagnetic Hamiltonian and study the dynamics of domain wall excita- 
tions in the presence of transverse magnetic field h x . We obtain dynamical spin 
correlation functions along the magnetic field S xx (q,uj) and perpendicular to it 
S yy (q,uj). It is shown that the line shapes of S xx (q,uj) and S yy (q,oj) are purely 
symmetric at the zone-boundary. It is observed in S yy (q,uj) for ir/2 < q < tt 
that the spectral weight moves toward low energy side with the increase of h x . 
This model is applicable to study the spin dynamics of CsCoCbj in the presence 
of weak interchain interactions. 



1 Introduction 

The spin-| Ising-like antiferromagnetic (AFM) chain has been the subject of theoretical 
studies for quite some time. The spin dynamics of the system are characterized by a 
picture of propagating domain walls or solitons. The magnetic compounds CSC0CI3 and 
CsCoBr3 are good examples of S — \ Ising-like AFM chains. The simplest exchange 
interaction Hamiltonian describing these compounds is the S = \ XXZ Heisenberg 
model 

#xxz = 2J]T [S*S z i+1 + e (S?S* +1 + SfSf +1 )} , < e < 1. (1) 

i 

For very small e, the lowest order ground state of Eq.(|l|) is the Neel states with a z- 
component of the total spin given by S , ^=0. Villain M has calculated the longitudinal 
correlation function S zz (q, uj) based on the basis states consisting of a single domain 
wall and predicted the appearance of a central peak with sharp shoulders. On the other 
hand, Ishimura and Shiba || proposed a picture of domain wall pair (DWP) states 
and showed that the propagating DWPs give rise to an excitation continuum around 
the Ising excitation energy 2 J. The transverse correlation function S xx (q, uj) exhibits a 
broad peak around 2 J. The existence of these peaks of S zz (q, uj) and S xx (q, uj) has been 
verified by inelastic neutron scattering experiments on CSC0CI3 || [|, ^ and CsCoBr 3 
0. A significant feature of the spin wave response of S xx (q, uj) near the zone centre 
(q = ir) is that the spectral weights are heavily concentrated towards the lower energy 
region. Nagler et. al. || added a staggered field term 

H s = hY / (-iyS? (2) 

i 

to the Hamiltonian in Eq.(jl]). The staggered field h has two contributions h Q and 
h lc . The first contribution originates from taking account of the exchange mixing of 
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higher levels with the ground doublet. The second contribution arises from the inter- 
chain exchange interactions at low temperatures. The interchain interactions treated 
in the mean-field approximation, give rise to the staggered field term h\ c . The effective 
Hamiltonian contains both the terms Hxxz and H$. With this effective Hamiltonian, 
the broad peak is found to spilt into discrete peaks which is known as Zeeman lad- 
der and observed in Raman scattering on CSC0CI3 and CsCoBr 3 J7j. However, the 
obsereved line shapes of S xx (q,u) are quite different from those of the theoretical pre- 
dictions. Matsubara el. al. || have included a weak next nearest neighbour (NNN) 
ferromagnetic (FM) interaction Hp in the Hamiltonian Hxxz in Eq. ([!]): 

Hf = -2J' £ [s*S? +2 + e (S*S? +2 + S y S y +2 )] . (3) 

i 

They have shown the existence of bound states of DWPs as well as the free DWP 
states and the transverse correlation function S xx (q,u) exhibits a sharp peak at lower 
energy region. The effect of transverse magnetic field h x on the spin dynamics of this 
model has been studied by Murao el. al. and shown that the spectral weight moves 
towards the low energy side in S yy (q, uj) for 7r/2 < q < n with the increase in h x , while 
there is no appreciable change in S xx (q,uj) for all q. The distribution of intensities of 
the sharp peaks in S yy (q,uj) vary irregularly for q w it. Although the proposed form 
of NNN FM coupling provides a good description of most of the experimental results, 



the required magnitudes of the NNN exchange \ J'\ ~ 0.1|J| is unphysically large [[10 



In 1996, Bose and Ghosh [Oj] have proposed the Ising-like fully anisotropic Heisenberg 
AFM Hamiltonian in ID and shown that the asymmetric line shapes of S xx (q,uj) and 
the bound states of DWPs can be derived. 

In the absence of the magnetic field, is a good quantum number and the eigen- 
values of different having unequal number of DWPs form different energy bands 
separating by energy 2 J. In the presence of longitudinal magnetic field, h z , still is a 
good quantum number and the eigenvalues within the same value of S% as well as the 
position of the peak of S xx (q, uj) shift parallel with the increase in h z . However, S% is 
no longer a good quantum number in presence of a transverse magnetic field h x and a 
mixing of states with different occurs. Thus, eigenvalues as well as eigenstates will 
be modified by h x and the characteristics of the spin dynamics will be different. 

In this paper, we study the effect of h x on the dynamical spin correlation functions 
in a fully anisotropic Ising-like S = \ Heisenberg AFM chain at low temperatures. 
Dynamical correlation functions S xx (q,uj) and S yy (q,uj) have been derived using the 
picture of propagating DWPs. Finally we introduce this model to explain the spin 
dynamics of CSC0CI3 taking into account the weak interchain interactions (Eq.(^)). In 
Section 2, the theory and the results for the eigenvalues of the DWP continuum and 
DWP bound states are derived. The dynamical spin correlation functions of CSC0CI3 
are presented in Section 3. Section 4 contains a discussion of the results obtained. 
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2 Model and domain wall pair states 



The one dimensional fully anisotropic Ising-like Heisenberg Hamiltonian in the presence 
of transverse magnetic field is given by 

N N 



h — 2^2 [JxSfSf +1 + j y sysf +1 + J Z S-S- +1 - sLL/iB-ffx E ^ 



N 



N 



2J ^2 S'i + 2 + 'S'j+i) + 2 + $i Si+l) E S'i 

i i 

J = J z , e x = Jx ^j Jy i e 2 = Jx ' 2J Jy ' , /i* = g±^BH x , e u e 2 < 1. (4) 

H x is the transverse magnetic field and assume h x <C 2J. N is the total number of 
spins. Since we are interested in excitations at low temperatures, we consider low 
lying excited states. These states can be obtained from the Neel state by flipping a 
block of adjacent spins, giving rise to DWP states with S^=0 and ±1 (Fig. 1). These 
excitations occur around the Ising energy 2J above the ground state. Following the 
method introduced by Murao et. al. 0, we classify these states into two series. Series 
'a' starts from the state with 5*^=1 where two domain walls are adjacent. Let m be 
the number (odd) of sites between two domain walls and (j^f 3 (m) be the corresponding 
Ising state. The subsequent states m a \m) (j = 2, 3, 4, • • •) are generated from 0^ (m) 
such that the separation between the domain walls is increased by unit lattice distance 
successively towards the right hand side of the chain. Hence, 



(5) 



| Neel) is one of the Neel states. We choose a linear combination of these basis states 
for describing propagating DWPs with wave vector q as 

b\?). = t/f E e-^f\m) (6) 

' m=odd 

On the other hand, series 'b' originates from the state with S% = —1, and the subse- 
quent states with S^=0 and —1 appear alternately, 

<f>?(m) = S-\Neel) S$ = -1, 

4>f{m) = S+y^fiH, (j = 2, 4, 6, • • •) ^ = 0, (7) 
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Taking linear combination of these states with wave vector q 

\j,Q)h = \hz E e-^f\m) (8) 



m=even 
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Figure 1: Neel states and DWP states for = ±1 and 0. The dotted vertical lines 
indicate the position of domain walls. 



With the help of the Eqs.(5)-(8), one can obtain H\j,q) a as: 

H\l,q) a = 2J\l,q) a + V ei \3,q) a + V e2 \l,q) b -^(\2,q) a + e^\2,q) b ) 



H\2, q) a = 2J\2, q) a + V ei |4, q) a - y q) a + |3, q) a + e~^|3, q) b + e^l, q), 



H\j,q) a = 2J\j,q) a + V ei \j+2,q) a + V*\j-2,q) a -^-(\j-l,q) a + \j + l,q) a ) 



(e*\j-l,q) b + e- i *\j + l,q) b ), j > 3, 



(9) 



where V ei = eiJ(l + e 2iq ) and V e2 = 2e 2 Jcosg. In the same manner, one could 
derive similar set of equations for H\j, q) b in terms of \n, q) a and \n, q) b . To avoid the 
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mixing between the states of series 'a' and 'b', we further introduce symmetric (a) and 
antisymmetric (/3) functions || defined as 

\j,q) a = -j=(\j,q)a+\j,Q)b), 

\U)p = -^(\j,q) a -\j,q)b). (10) 

Hence, one can express H\j,q) a as 

H\l,q) a = (2J + V e2 )\l,q) a + V ei \3,q) a -V a \2,q) a , 
H\2,q) a = 2J\2,q) a + V ei \A,q) a -(V:\l,q) a + V a \3,q) a ), 

H\j,q) a = 2J\j,q) a + V ei \j+2,q) a + V e * 1 \j-2,q) a -(V:\j-l,q) a + V a \j+l,q) a ) 

J>3, (11) 

where V a = 4f (1 + e~ 4<? ). Similarly, one can derive H\j,q)p with a and V a being 
replaced by (5 and Vp — (1 — e~ iq ), respectively The first excited states can be con- 
structed as a linear combination of symmetric and antisymmetric functions separately, 

*«(g) = X>ib'>?>a and ^ p {q) = ]T/? J |j,g) /3 . (12) 
i i 

With the help of the Eq.([ll|), the following equations for the coefficients acj and f3j are 
obtained as 

\ a ax = (2J + V e2 )ai + 140:2 + V ei a 3 , 
X a a 2 = 2Ja 2 + V a («i + a 3 ) + V^a^, 



A a «j = 2J«j + V a (oj-x + Oj+i) + V r ei (aj_2 + 0^+2) > J > 3, (13) 

where A a is the eigenvalue, V" = —h x cos (|Y Ki = 2eiJcosg and «j = In the 

same manner, one can derive similar set of equations for (5 with a being replaced by 

( \ _ — i(fc + 7r) . 

|J and ctj by f3j = fye a 

Dispersion relations are obtained numerically by solving Eqs.fllB"!), with N=1000. 
Here, we present the results for e\ = 0.05 and t 2 = 0.10, since these values are esti- 



mated in the compound CSC0CI3 ||11|| . Figure 2 shows the disperson relations in the 
symmetric and antisymmetric modes, X a and Xp, respectively, for h x = 0. The spin 
wave continuum and the bound state energy are plotted by solid and dotted lines, 
respectively. In the symmetric mode, the bound state energy lies above the continuum 
for q< \ and below the continuum for q > | , while the reverse is true for the antisym- 
metric mode. The bound state does not exist when €\ > e 2 . When h x 7^ 0, the energy 
band extends towards the high energy region for < q < | in both the symmetric and 
the antisymmetric modes. The spin wave excitations have a width at q = | in contrast 
with the case of h x = 0. The width also broadens with the increase of h x . Note that 
the bound state energy is not affected by the presence of h x . 
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Figure 2: Spin- wave excitation continuum (solid lines) and DWP bound-state energies 
(dotted lines) of the symmetric and antisymmetric modes for ei = 0.05 and e2 = 0.1. 



3 Dynamical spin correlation functions at T=0 K 

The dynamical spin correlation functions along the direction of h x at T=0 is defined 
S"(q,u,) = El(*c \S x (q)\ V g )\ 2 6(uj - A c + A g ) (14) 



as 



where |^ g ), |^ e ) denote the ground and excited states, respectively, and A g , A e are the 
corresponding eigenvalues. In this case, the ground state is one of the Neel states and 
the summation extends over the first excited states only Also, the Fourier transform 
of S*, i.e., 



S x (q) = V e 4 ^ (Sf + S- 



2VN , 

Similarly, the dynamical spin correlation function perpendicular to the direction of h x , 
S yy (q,uj), is defined by replacing the superscript x with y in Eq. (|Tjj) , where 

S y {q) = V e 1 ^ (Sf - Sr 

Since the ground state is the Neel state, S xx (q, uj) and S yy (q, ui) directly reflect the wave 
number dependence of the excited states. With the help of the Eq.(^), the dynamical 




Figure 3: The functions S xx (q,u) and S yy (q,cu) for q — 0, f and 7r in the histogram 
with Aco = 0.025J. 



7 



spin correlation functions can further be written as 



= j£K| 2 ^-A Q + A g ), 

S^faw) = - £ |ft| 2 *(a; - + A g ). (15) 

Note that S xx (q, to) depends only on |«i| 2 while S yy (q, to) on |/9i| 2 . Thus, the symmetric 
mode is directly reflected on S xx (q, to), whereas the antisymmetric mode on S yy (q,io). 
The functions S xx (q,to) and S yy (q,to) for h x = are shown in figure 3. A sharp peak 
originates from the bound state, whareas the broad peak from the free DWP states. 
The intensity of the sharp peak does not depend on the number of spins N, while the 
broad peak comprises (N — 1) peaks of which has intensity of the order of 4. Note 
that at the zone-boundary (q = |), the width of the continuum vanishes. This is also 



verified in neutron scattering experiments on CSC0CI3 [|12 |. 

The line shapes of S xx (q, to) at h x = 0.1J have been plotted in figure 4. The main 
feature of S xx (q, to) induced by h x shows that the line shape is purely symmetric at the 
zone boundary and it is highly asymmetric away from the zone boundary. At q ~ 0, the 
sharp peak occurs at high energy and the tail at lower energy region. For < q < |, 
the spectral weight concentrates mainly in the middle of the continuum. The sharp 
peak emerges again in the lower energy region for | < q < tt. At q « 7r, the line shape 
is not affected by h x as expected from the dispersion relation shown in figure 2. Figure 
5 shows S yy (q,io) at h x = 0.1J. The line shapes of S yy (q,io) is again symmetric at 
q = ~. The q dependence of S yy (q, 10) is opposite from that of S xx (q, 10) and it remains 
unaffected at q = 0. For < q < ~, the sharp peak appears on the lower energy 
region of the broad peak. The sharp peak originates from the bound state, whereas 
the broad peak from the DWP continuum. For | < q < tt, the sharp peak appears on 
the higher energy side of the broad peak and the tail is found to enhance towards the 
lower energy region. In figures 6 and 7, S xx (q, to) and S yy {q,uj) are shown for different 
values of h x , respectively. With the increase of h x , the features mentioned above are 
enhanced. The height of the sharp peak is found to diminish with the increase of h x . 
Note that S xx (q, 10) for < q < | and S yy (q,io) for | < q < n are sensitive on h x as 
observed by Murao et.al. M. 



3.1 Interchain interactions in CSC0CI3 

Now we add the staggered field term to the Hamiltonian in Eq. (^) . The full Hamiltonian 
looks like 



N T ei e 2 



N 



N 



- hx j2s x -h ic j2(-iys?. 



(16) 



Here, we further consider that the staggered field h ic originates due to the weak in- 
terchain interaction. The interchain interaction has been treated in the mean-field 
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Figure 4: The function S xx (q, uj) for different values of q and h x = 0.1J. The width of 
the histogram is Atu = 0.025J. 
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Figure 5: The function S yy (q,ui) for different values of q and h x = 0.1J. The width of 
the histogram is Au; = 0.025J. 
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Figure 6: The function S xx (q, uj) for different values of q and h x . The width of the 
histogram is Auj = 0.025J. 
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12 



approximation. Thus, 

h ic = 2J'Y,(S? +5 ), (17) 

<5 

where J' is the interchain interaction strength and S is the nearest neighbour on the ab 
plane. Following the same technique developed in Section 2, we obtain S xx (q,uj) and 
S yy (q,uj) numerically. 

In figure 8, S xx (q, uj) for different values of h- lc are shown. When h ic ^ 0, the broad 
peak of S xx (q,u) has been splitted into discrete peaks, which are known as Zeeman 
ladder jl3| . With the increase of h ic , the number of peaks decreases and the separation 
between them becomes wider. For | < q < ir, the intensity of the peak at the lowest 
energy is stronger. Several peaks having nearly the same intensity are observed at 
q « 0. Figure 9 shows S yy (q,u) for different values of h- lc . For < q < |, several 
strong peaks appear in the lower energy region. Note that the difference between 
S xx (q,uj) and S yy (q,u) are remarkable at q = n. In S yy (q,uj) for q = n, several strong 
peaks appear in the higher energy region and weak peaks in the lower energy region, 
and their intensities vary irregularly with the increase in h lc . Thus the combined effect 
of h x and the interchain interactions is also found in S xx (q, uj) and S yy (q,uj) for q > | 
as reported in the Ref. ||. 

In CSC0CI3, the magnetic Co +2 ions are surrounded by trigonally distorted octahe- 
dra of CT~ ions and form chains along the c-axis with successive octahedra sharing a 
common face. The CsCl^ chains are arranged in a triangular array. Since the exchange 
coupling between chains is antiferromagnetic, the triangular array forms a frustrated 
system. Thus there is no possibility of a perfectly regular AFM ordered state. Two 
different three-dimensionally ordered phases occur in CSC0CI3. First, below ~ 
21K partially disordered AFM phase (A) is formed in which one-third of the chains 
are paramagnetic. A phase change takes place below T/v 2 ~ 10-14K, to a ferrimagnetic 
phase (F), in which the paramagnetic chains align in the same direction, so that two- 
thirds of the chains are aligned in one direction and one-third in the opposite direction 
(Fig. 10) 0. 

It is obvious from the Eq. (|TT|) , that the staggered field h ic could take one of the 
two possible values of h lc = 6 J' and 0. Similarly, in the partially disordered phase, 
h ic takes one of the four values of h ic = 6 J', 4J', 2 J' and 0. The functions S xx (q, uj) 
and S yy (q,uj) of CsCoCl 3 are obtained by summing up those of the individual chains. 
In the ferrimagnetic phase, ratio of the number of chains with h ic = QJ' and are 
estimated as 1 : 2. In the partially disordered phase, that with h ic = 6J', 4J', 2J' 
and are estimated as 1 : 3 : 3 : 5 0. Figure 11 shows S xx (q, uj) for q = and 
S yy (q,uj) for q = it in the ferrimagnetic phase. The line shape of S xx (q,u) is affected 
strongly by the h x . As h x increases, the difference in intensities among various discrete 
peaks reduce markedly and the sharp peak disappears. On the other hand, S yy (q, uj) is 
affected slightly by h x . The spectral weight enhances towards the lower energy region 
with the increase in h x , and the intensity distribution varies irregularly. S xx (q, uj) for 
q = and S yy (q,uj) for q = it obtained in the partially disordered phase are shown in 
figure 12. S xx (q,uj) is strongly affected by h x while S yy (q,uj) is much less. S xx (q, uj) 
for q = tt and S yy (q,uj) for q = are not affected by h x both in the ferrimagnetic and 
the partially disordered phases. 
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Figure 8: The function S xx (q, co) for different values of h K and h x = 0.1J. The width 
of the histogram is Aco = 0.025J. 
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Figure 9: The function S yy (q,cu) for different values of h ic and h x = 0.1 J. The width 
of the histogram is Auj = 0.025J. 
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(a) (b) 



Figure 10: Magnetic ordering in ab plane below Tjy^: (a) Ferrimagnetic structure for 
T < T Na (F-phase). (b) Structure for T Na < T < T Nl (A-phase). Chains marked with 
open circles are disordered. 

4 Discussion of results 

We have studied the effect of the transverse magnetic field h x on dynamical properties 
of one dimensional fully anisotropic Ising-like antiferromagnet at low temperatures. We 
have shown using this Hamiltonian that some of the results obtained by Murao el. al. 
0, in which a FM NNN interaction is assumed besides the usual AFM NN interaction, 
can be qualitatively reproduced. These include the formation of DWP bound states, 
two types of excited modes which are symmetric and antisymmetric with respect to 
the states with = 1 and —1 and an asymmetry in the line shape of the correlation 
functions S xx (q, u) and S yy (q, to). In order to obtain the asymmetry in the line shape of 
S xx (q,ui) and S yy (q,uj), if a FM NNN exchange of magnitude \ J'\ ~ 0.1|J| is required, 
which considering that the NNN exchange is through two nonmagnetic ligands would 
seem to be unphysically large |ID) . On the other hand, our model could explain all these 
charateristics with the usual NN AFM exchange interactions. There are, however, a 
number of differences. Murao et. al. observed a single bound state branch which is 
symmetric with respect to the zone boundary, whereas in the present study, different 
bound state branches are obtained for symmetric and antisymmetric modes which are 
asymmetric with respect to the zone boundary. No experimental evidence are as yet 
available on the effect of bound states on the thermodynamic and dynamic properties of 
the compounds CSC0CI3 and CsCoBr 3 . The symmetric modes contribute to S xx (q,u>), 
whereas antisymmetric modes contribute to S yy (q,u). Both S xx (q,u) and S yy {q,uj) 
have symmetry at the zone boundary even in the presence of either both h x and h lc 
or any one of them. This symmetry is totally lost away from the zone boundary. This 
theory is also valid to analysis of spin dynamics of CSC0CI3 and CsCoBr3 at finite 
temperatures, because the correlation length of the spin along the chain is very large 
even in the paramagnetic phase (T « 21K). Thus the effects of h x on S xx (q,u) and 
S vy (q,ui) discussed here are yet to be observed in a real system. 

Apart from relevance to experimental systems such as CsCoCl 3 and CsCoBr 3 , 
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Figure 11: The function S xx (q, uj) for q = and S yy (q,uj) for g = 7r in the ferrimagnetic 
phase. The width of the histogram is Au; = 0.025J. 
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Figure 12: The function S xx (q, oj) for q = and S yy (q,oj) for q = n in the partially 
disordered phase. The width of the histogram is Aoj = 0.025J. 
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the present study is intended to provide insights about the spin dynamics of fully 
anisotropic Ising-like AFM system in the presence of transverse magnetic field h x . 
The ground state energy and low lying excitation spectrum of the fully anisotropic 
Hamiltonian are known exactly because of the mapping between the fully anisotropic 



Hamiltonain and the exactly solvable eight vertex model [16]. Our calculations pro- 
vide us with some physical insights about spin dynamics in Ising-like fully anisotropic 
AFM system in the presence of transverse magnetic field. 
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